Alternative mechanisms of structuring biomembranes: Self-assembly vs. self-organization 



O 

o 

(N 
> 

o 



U 
d 

> ■ 

m 
cn : 
o ■ 
\o ■ 
o 
vn 
o 



I 



X 



Karin John' B and Markus Bar^'0 
' Max-Planck-Institut fiir Physik komplexer Systeme, Nothnitzer Str. 38, D-01 187 Dresden, Germany 
^ Physikalisch-Technische Bundesanstalt, Abbestr. 2-12, D-10587 Berlin, Germany 

We study two mechanisms for the formation of protein patterns near membranes of living cells by mathemat- 
ical modelling. Self-assembly of protein domains by electrostatic lipid-protein interactions is contrasted with 
self-organization due to a nonequilibrium biochemical reaction cycle of proteins near the membrane. While 
both processes lead eventually to quite similar patterns, their evolution occurs on very different length and 
time scales. Self-assembly produces periodic protein patterns on a spatial scale below 0. 1 [jm in a few seconds 
followed by extremely slow coarsening, whereas self-organization results in a pattern wavelength comparable 
to the typical cell size of 100 |jm within a few minutes suggesting different biological functions for the two 
processes. 

PACS numbers: 87.18.Hf 87.15.Kg 82.40.Ck 64.75.+g 



Living cells display internal structures on various lenght 
scales, that are regulated dynamically. Examples include the 
organisation of nerve cells into axon, body and dendrites iQl] 
or the occurrence of short-lived signaling patches in the mem- 
brane of chemotaxing amoebae The origin of many of 
those structures is a nonuniform distribution of biochemical 
molecules, that can be achieved by a spontaneous symmetry 
breaking through local fluctuations. The diversity of time and 
length scales in cellular structures strongly suggest that a vari- 
ety of mechanisms participate in the structuring process. From 
a physics perspective, one can distinguish at least two differ- 
ent classes of processes: self-assembly and self-organization 

a. 

Self-assembly implies spatial structuring as a result of min- 
imization of the free energy in a closed system. Hence, a self- 
assembled structure corresponds to a thermodynamic equilib- 
rium. A prominent example for self-assembly in single cells is 
phase separation of lipids and proteins due to macromolecular 
interactions. Phase separation occurs if the interaction ener- 
gies dominate the entropy contribution. Model membranes 
show a phase separation due to lipid-lipid interactions |0, [Hi 
or lipid-protein interactions [6j] (for a critical discussion see 
also [7]). Theoretical analysis of phase separation in biomem- 
branes has been mainly restricted to free energy considera- 
tions |8, 9], which can predict the final equilibrium state, but 
do not capture the transient dynamics. 

Self-organization, in contrast, requires a situation far away 
from thermodynamic equilibrium and is possible only in open 
systems with an external energy source. One prominent exam- 
ple is pattern formation in reaction-diffusion systems, which 
has first been proposed by Turing [ 10^ and later on become in- 
fluential in development biology 1 11]. The key ingredients are 
nonlinear self-enhancing reactions, a supply of chemical en- 
ergy and competing diffusion of the involved molecules. Ex- 
perimental patterns in single cells that are successfully mod- 
elled by reaction-diffusion equations include calcium waves 
lfl2ll and protein distributions in E. coli li 1 3|] . 



In the vicinity of membranes molecular interactions and 
reaction-diffusion processes occur simultaneously. In this 
Letter we discuss and compare two dynamical models for 
the alternative mechanisms of pattern formation near the 
cell membrane: model I describes self-assembly of protein 
domains due to lipid-protein interactions and model II de- 
scribes an active reaction-diffusion mechanism resulting in 
self-organization of proteins. 

The two models are based on the properties of the GMC 
proteins reviewed in lfl4l[l5l] . which play a critical role in the 
development of the neural system [16] and in the regulation 
of cortical actin-based structures and cell motility |17]. A 
key result is that both models display similar qualitative dy- 
namical behavior but act on largely different time and length 
scales. Model I leads to the formation of stationary domains 
on a submicrometer scale within seconds. A slow coarsen- 
ing process follows the expected power law. Phase separa- 
tion arises from a relaxational process into thermodynamic 
equilibrium, where the equilibrium state is characterized by 
vanishing fluxes. Model II leads to the formation of large 
stationary structures on the scale of a eukaryotic cell within 
lOmin. The steady state of this active phase separation (= 
phase separation from self-organization) is characterized by 
nonzero fluxes. 

Model I is based on attractive interactions of GMC pro- 
teins with acidic lipids in membranes 1 18]. In this mechanism 
(Fig.[T](a)) proteins (area fraction Cm) are associated with a 
lipid membrane, consisting of type 1 (area fraction c/) and 
type 2 (area fraction 1 — c/) lipids of identical size. The size 
ratio between proteins and lipids is N . Proteins interact attrac- 
tively with type 1 lipids. Proteins and lipids are allowed to dif- 
fuse in the plane of the membrane. The number of bound pro- 
teins and the average membrane composition are conserved 
quantities. We write the free energy of the protein covered 
membrane similar to ]8] as 



F = keTNa / dA 



(1) 
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where dA denotes integration over the membrane surface, Na 
denotes the number of lipids per unit area and 

// = Q InQ + (1 - Q)ln(l - q) + 
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The local part of the free energy ^ consists of the entropic 
contributions of the hpid (first and second term) and the pro- 
tein phase (third and fourth term) and the interaction energy 
between lipids and proteins (last term) [19]. For simplicity 
we assume that electrostatic repulsion and non-electrostatic 
attraction in the lipid and protein phase cancel locally. This 
is somewhat arbitrary, but does not change the nature of the 
results as long as demixing is governed by lipid-protein inter- 
actions. The non-local part of the free energy ([T]) is governed 
by the interfacial energy in the lipid phase, which is a small 
quantity. For the parameter x we set x ~ \uiP from Cahn- 
Hilliard Theory [20] where ui and I are the typical interaction 
energy (on the order of several fee 7") and length (on the or- 
der of a lipid headgroup size, i.e. / « 1 — 2 nm), respectively. 
Evolution equations for c; and Cm are obtained using linear 
nonequilibrium thermodynamics and the mass balance equa- 
tions to give 



dtci = -V • ji with ii = 

dtCm = -V • jm with im = 



(3) 



More specifically, the fluxes ji and jm are given by 

J, = -Di [Vci - q(1 - q) (uVc,„ + xV^q)] (4) 

jm [Vc„ - NCm{l " C„)uVCi] , (5) 

where Di = MikeTNa and Dm - MmksTNa. 
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FIG. 1: Schematic depiction of models I (a) and II (b). (a) The mem- 
brane consists of a mixture of type 1 (black) and 2 (white) lipids with 
adsorbed proteins, (b) Proteins undergo a cycle of phosphorylation 
and dephosphoylation as explained in the text. 



Model II is based on a biochemical cycle of GMC phospho- 
rylation and dephosphorylation, called myristoyl-electrostatic 
(ME) switch 1I2III . In contrast to model I, we will assume 
that all different lipid species are distributed uniformely in 
the membrane. In this phase separation model (Fig.[T](b)) a 
protein can associate reversibly with the membrane following 
a mass action law. On the membrane proteins are irreversibly 
phosphorylated by a protein kinase, which disrupts membrane 



binding immediately and translocates the protein into the cy- 
tosol where it is dephosphorylated by a phosphatase. We can 
cast these processes into a three-variable reaction-diffusion 
model of the following form 



dtCm — ^ad(l (^m)Cc ^rfeCr? 

fcfci(l ^ Cm)-, ; 



Dm^ Cn 



(6) 



df Cc — kdeCm ~ kad{^ Cm)Cc + kphCp + Z?cV^Cc (7) 



dtCp 



^fci ( 1 Cm ) ~ 



kphCp + Dc'^'^Cp . (8) 



Cm, Cc and Cp denote the concentrations of membrane bound, 
cytosolic unphosphorylated and cytosolic phosphorylated pro- 
teins, respectively, had and kde denote the rate constants of 
membrane association and dissociation of the unphosphory- 
lated protein and k^i and kph denote the enzymatic activities 
of the kinase and phosphatase. For the kinase activity we have 
used a Michaelis-Menten type kinetics ll22ll . whereas we have 
neglected this property for the phosphatase, assuming that the 
concentration of phosphorylated proteins is well below the re- 
spective Michaelis-Menten constant |23]. Based on the prop- 
erties of protein kinase C we assume, that the kinase needs 
lipids for full activation ll24l I25I1 and that membrane bound 
proteins decrease the available membrane space and thus the 
kinase activity [6]. Eqs. (|6]l-([8j constitute a non-equilibrium 
system, since the kinase activity is sustained and consumes 
ATP which is produced by the metabolism. The total protein 
concentration ct — y J dV {cm + Cp + Cc) is conserved. 

The parameters for both models have been taken from ex- 
periments. First we consider the linear stability of the uniform 
steady states (Fig.|2]i. In both models we can identify a region 
of linear instability, characterized by real positive eigenvalues. 
Using the Maxwell construction and tie-lines we calculated 
for model I the equilibrium state for a given uniform state and 
identified a metastable region, where a stable uniform state 
coexists with a stable demixed state. The demixed state is 
characterized by a phase with a high concentration in protein 
and lipid 1 and a phase with a low concentration in protein and 
lipid 1 . Using continuation methods we computed stationary 
solutions for model II. Fixing the wavelength to 120 |am (com- 
parable to the size of a eukaryotic cell) one can also identify 
regions, where a linearly stable uniform solution coexists with 
a linearly stable stationary periodic solution. 

Figs.[3](a) and (b) show the largest eigenvalues for both 
models from the linearly unstable regions in Figs.|2](a) and 
(b). The instabilities belong to the type lis class l29tl . which is 
characterized by a real critical eigenvalue with wave number 
zero and can be attributed to the conservation relations in both 
models. Although both models have instabilities of the same 
type they develop on very different length and time scales. 
The wavelength of the fastest growing mode Am in model I is 
linked to the molecular interaction length. In our example in 
Fig.[3]a it is on the scale of 50 nm with a growth rate of 10 s^^. 
In contrast, A,„ of model II is determined by kinetic rate and 
diffusion constants and is of the order of 10 |jm. In the exam- 
ple in Fig.[3]b the corresponding growth rate is 0.01 s^^. 

The results of the linear stability were confirmed by nu- 
merical simulations in two dimensions with periodic bound- 
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FIG. 2: (a) Phase diagrams with tie-lines for the phase separation 
in model I and (b) phase diagram for the active phase separation in 
model II. Parameters in (a) are A'^ = 25 and u — 1.0 and in (b) are 
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FIG. 3: Dispersion relations for models I (a) and II (b) for the uni- 
form solution. Shown is only the largest (real) eigenvalue. Pa- 
rameters in (a) are A'^ — 25, C; = 0.2, Cm = 0.5, u = 1.01, 



Di = 1.0 /xm-" 
and in (b) Ct = 
in Fig.|l](b). 
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ary conditions. Simulations were started from the uniform 
steady state with small amplitude perturbations. Figs.|4](a) 
and |4](b) show exemplary two simulations for models I and 
II with parameter values from the linearly unstable regions 
in Fig|2] In both models stationary structures develop, which 
are not stable but display a coarsening behavior for later 
stages. For model I we found the scaling law k{t) — at^^ 
(Fig-EJ with an exponent v k, 1/A consistent with a mod- 
ified Lifshitz-Slyosov- Wagner theory [31] for concentration 
dependent mobility coefficients. Since we are considering the 
two-dimensional case one should rather expect a growth law 
^(t) = ar^" with a modified time scale r = t/ ln(4t) ll32ll . 
However, on the time scale of our numerical experiments both 
growth laws yielded similar exponents. Although the initial 
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FIG. 4: Snapshots of a two-dimensional numerical simulations of 
models I (a) and II (b) with periodic boundary conditions. Shown is 
the concentration of membrane bound protein Cm . Parameters in (a) 
and (b) are as in Fig.[3](a) and (b), respectively. 
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FIG. 5: Late stage behavior of model I (t* = tDi/x, k* — k,^/x)- 
Shown are loglog-plots of the evolution of the mean wave number 
of the structure factor of the membrane bound protein concentration 
with time (a) or a modified time (b) as described in the text. The 
exponent v of the power law was obtained by a least square fit. Pa- 
rameters are the same as in Fig.[3](a). 



growth is very fast the initally observed wave numbers are 
small. To obtain structures on the scale of the cell coarsening 
has to occur over several orders of magnitude with v = 1/4, 
which is a slow process. The possible scaling behavior of 
model II is irrelevant for practical purposes since the initially 
developed structures are already on a scale comparable to the 
system size and the first or second coarsening step will lead to 
polarized cells with a single domain of high concentration of 
membrane bound proteins. One can easily see in Fig.|4](b) that 
large structures have appeared after ten minutes and within 
one hour coarsening is complete. 

In this Letter we have introduced and analyzed two alterna- 
tive models for pattern formation of GMC proteins. GMC pro- 
teins are on the one hand found to form domains by virtue of 
their attractive electrostatic interaction with acidic lipids from 
self-assembly. On the other hand, they can exploit an ATP- 
driven phosphorylation-dephosphorylation cycle (myristoyl- 
electrostatic switch) for their self-organization. The striking 
difference between both mechanisms lies in the relevant time 
and length scales. The spinodal length scale of the phase sepa- 
ration in model I is closely linked to the molecular interaction 
length in the lipid phase. For fluid membranes under physio- 
logical conditions the relevant interaction scale is comparable 



4 



to the size of a lipid molecule (« 1 nm). The initial struc- 
ture formation is fast with growth rates of the order of 10 s^^, 
but the coarsening process follows a scaling law. Thus coars- 
ening does not lead to structures on the size of the cell in a 
biologically relevant time. The reaction-diffusion mechanism 
(model II) leads initially to large structures, which are on the 
scale of an eukaryotic cell. A rough estimate for the length 
scale is given by the quantity V Dt, where D w lOpm^ s^^ 
is the typical intracellular diffusion constant of a protein and 
T = 10 s is a typical time for a biochemical reaction. This 
yields a length scale of « 10|jm. A combination of both 
mechanism leads to a more complicated model, that displays 
oscillatory dynamics and traveling domains 13311 . Here, we 



have used specific physical and chemical properties of GMC 
proteins, but the typical scales for molecular interaction en- 
ergies and ranges as well as for reaction rates and diffusion 
constants will be comparable for other processes near mem- 
branes. We propose that both mechanisms are relevant for 
different aspects of structuring membranes: protein-lipid in- 
teractions are suitable for rapid structuring of membranes on 
a submicrometer scale, whereas reaction and diffusion of pro- 
teins produce a structure on the scale of the size of typical 
eukaryotic cell within minutes and are potentially useful for 
polarizing a whole cell into two main compartments (e. g. 
front and back). 
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